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We study the cosmic microwave background temperature and polarisation spectra sourced by 
multi-tension cosmic superstring networks. First we obtain solutions for the characteristic length 
scales and velocities associated with the evolution of a network of F-D strings, allowing for the 
formation of junctions between strings of different tensions. We find two distinct regimes describing 
the resulting scaling distributions for the relative densities of the different types of strings, depending 
on the magnitude of the fundamental string coupling gs . In one of them, corresponding to the value 
of the coupling being of order unity, the network's stress-energy power spectrum is dominated by 
populous light F and D strings, while the other regime, at smaller values of g^, has the spectrum 
dominated by rare heavy D strings. These regimes are seen in the CMB anisotropies associated 
with the network. We focus on the dependence of the shape of the B-mode polarisation spectrum 
on gs and show that measuring the peak position of the B-mode spectrum can point to a particular 
value of the string coupling. Finally, we assess how this result, along with pulsar bounds on the 
production of gravitational waves from strings, can be used to constrain a combination of gs and the 
fundamental string tension ^f- Since CMB and pulsar bounds constrain different combinations of 
the string tensions and densities, they result in distinct shapes of bounding contours in the {^F,gs) 
parameter plane, thus providing complementary constraints on the properties of cosmic superstrings. 

I. INTRODUCTION 

Although it has been established that a network of cosmic strings cannot source the majority of the observed cosmic 
microwave background (CMB) temperature anisotropy [1], the CMB can still provide a distinctive signature of their 
presence through the specific primordial B-modc polarisation spectrum [2-9]. The spectrum generated by strings is 
different from the one generically produced from tensor modes arising in inflationary scenarios, and future probes of 
the B-mode should be able to reveal the presence of cosmic strings, even if strings contribute as little as 0.1% to the 
CMB temperature anisotropy [5-9]. 

Interest in cosmic strings has revived following the realisation that they can arise in superstring theory [10, 11], 
for example in models of branc inflation [12-15] (for a review see [16]). Cosmic superstrings can have small tensions 
(10~^^>G'^>10~^ [10, 11, 17]), can be effectively stable over cosmological timescales, and can stretch over cosmological 
distances [18-20]. Hence, they can have interesting cosmological implications. Furthermore, their intcrcommutation 
probabilities can be significantly less than unity [11, 21-23] and, because of the charges present on them, they can zip 
together to form Y-junctions (trilinear vertices) , leading to more complicated networks than those usually considered 
in the case of 'standard' Abelian cosmic strings. Understanding the imprint of such additional network features on 
obscrvables, such as CMB temperature and polarisation, is a step that may lead to interesting new constraints on the 
basic parameters of string theory, such as the string coupling gs and the fundamental string tension ^p. 

Several approaches have been developed to model the evolution of cosmic string networks, and an interesting 
recent attempt to extend them to cosmic superstring networks - which contain different types of string - is due to 
Tye, Wasserman and Wyman [24]. Their model, based on the "velocity dependent one-scale" model of Martins and 
Shellard [25, 26], describes evolution of a multiple tension string network (MTSN) under the assumption that all types 
of strings have the same correlation length and root-mean-square velocity. By studying the evolution of the number 
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density of strings, they find that scahng is achieved when the energy associated to the formation of junctions is 
assumed to be radiated away. This model has been extended in [27] , where the authors assigned a different correlation 
length and velocity to each string type, and enforced energy conservation at each junction. Scaling is again achieved 
(with different number densities), but not as generically as in [24]. 

In a complementary approach, a number of authors have studied the kinematics of cosmic string collisions [28, 29]. 
When two Nambu-Goto (NG) strings (of generally different tensions) collide, rather than intercommuting in the 
standard way, they can form two junctions and a linking string of a third tension. Kinematically this can only occur 
if the relative orientation, velocity and string tensions lie in certain ranges. In [30], the authors extended their earlier 
studies to {p, g)-cosmic superstrings by modifying the NG equations to take into account the additional requirements 
of flux conservation. Once again the kinematic conditions required for the formation of Y-junctions were established, 
with results very similar to the ones obtained for NG strings. These kinematic constraints have been checked quite 
extensively with dynamical field theory simulations of strings collisions, and the agreement is (generally) good [31-35]. 
Recently, in [36] these constraints have been incorporated into the generalised MTSN velocity one-scale model of [27], 
leading to new conditions required for scaling and thereby providing the most complete model of cosmic superstring 
evolution to date. 

In this paper, we use the model of [36] to study the evolution of a cosmic superstring network for different values 
of the string coupling gs and different charges (p, q) on the strings. We find that in all cases the three lightest 
strings, i. e. the (1,0), (0, 1) and (1, 1) strings, dominate the string numfter density. When the string coupling is large, 
gs ^ 0{1), most of the network energy density is in the lightest (1, 0) and (0, 1) strings (respectively F and D strings), 
whose tensions are approximately equal. At smaller values of gs ^ 0(10"^), the (1,0) string becomes much lighter 
than both the (0,1) and (1,1) strings, and dominates the string number density. However, because of their much 
larger tension, the energy density of the network at small couplings can be dominated by the rarer (0, 1) and (1,1) 
strings. The existence of these two distinct limiting scaling behaviours at large or small values of gs is quite generic, 
although the specific details are somewhat dependent on the model-dependent value of the effective volume of the 
compactified dimensions. In either of the two limiting regimes, the energy density of the multi-tension network is 
effectively dominated by strings of one tension. 

With the scaling solutions to hand we then focus on the CMB imprints of these networks, using a modified version of 
the publicly available code CMBACT [38, 39]. In particular, we extend the Unconnected Segment Model (USM), first 
introduced in [1, 40], to describe the MTSN of [36] and implement it in CMBACT to obtain the CMB temperature 
and polarisation spectra. We find that for sufficiently large values of the parameter which is inversely proportional 
to the effective volume of the compactified dimensions (see Eq. (17) below), the two limiting regimes, one with the 
network energy dominated by light populous strings and the second with it dominated by rare heavy strings, can each 
produce distinct shapes for their CMB B-mode polarisation. In particular, for w ^ 1, the position of the peak in the 
B-mode spectrum is at £ w 770 for gs = 0.9 and at ^ sa 610 for gs = 0.04. This opens up the exciting possibility that 
upcoming observations may not only constrain the overall contribution of strings, but in fact rule out certain values 
of the string coupling. Namely, the combination of the normalisation and the peak position of the B-mode spectrum 
can point to a particular combination of gs and the fundamental string tension ^p- 

It is common to report constraints on standard cosmic strings in terms of bounds on the single dimensionless string 
tension G/i. These bounds have an implicit assumption on the number density of strings corresponding to the usual 
Abelian Higgs model strings with intcrcommutation probability V = I. However, in a more general situation such as 
that of cosmic superstrings considered here, there can be smaller intcrcommutation probabilities as well as strings of 
different tensions. As a result, each type of string will in principle have a different number density: the same fraction 
of CMB anisotropy can be sourced either with many light strings or with a few heavy ones. In general, each type 
of observational bound will constrain a different combination of the string tensions and densities (which, for cosmic 
superstrings, are derived from the fundamental string tension fip as well as gs). In particular CMB and pulsar bounds, 
which we discuss in Section IV, will lead to different shapes of bounding contours in the (fip^gs) parameter plane. We 
show that combining these two constraints can lead to complementary constraints on properties of superstrings. The 
position of the peak in the B-mode spectrum can be used to further eliminate a large region of the {fip, gs) parameter 
space. 

This paper is organised as follows. In Section II we summarise the extended VOS model which describes multi- 
tension networks with junctions, and for the case of cosmic superstrings we present the scaling solutions as a function 
of gs ■ In Section III we determine the temperature and B-mode spectra for these scaling solutions using a generalised 
version of CMBACT. In Section IV pulsar constraints on gravitational waves from string networks are discussed, and 
we conclude in Section V. 
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II. SCALING OF COSMIC SUPERSTRING NETWORKS 



Cosmic superstrings provide an example of a network of strings witli multiple tensions, and one in which the strings 
can join at Y-shaped junctions. A general (p, q) cosmic-string has p quanta of F charge (associated with fundamental 
strings) and q quanta of D charge (carried by D-branes) [17, 41, 42]. 

We begin this section by briefly reviewing the velocity dependent one-scale (VOS) model for single-tension strings, 
before generalising it (section II B) to a model for a network of TV different types of strings with junctions following 
[27, 36]. We then customise the parameters of this model to the case of {p,q) superstrings (section II C) paying 
particular attention to their dependence on the string coupling gs- While the exact values of the parameters are 
model-dependent, e.g. depend on the choice of the compactification manifold, we are able to identify general trends in 
their dependence on Qs. These trends, in turn, lead to two distinct scaling scenarios in the limits of large and small gs 
(section II D) that may be distinguished observationally. To gain additional intuition, we have included an Appendix 
where we analytically solve for the scaling parameters in the limit of small gs to complement the exact numerical 
solutions. 



A. The VOS model for a network of a single type of string with no junctions 

The so-called "one-scale model" [43] assumes that the string network can be characterised by a single length scale, 
namely the correlation length L defined by 



5. (1) 



where p is the energy density in the infinite string network and /i is the string tension. In principle, there are at least 
two different fundamental length scales in the network — the typical "smoothness" length L of long strings, and the 
average distance between strings L. The one-scale model takes the two lengths to be equal {L=L), an approximation 
which appears to be reasonably well satisfied [44] for a network of NG strings of tension /i with intercommutation 
probability V — I. 

As the network evolves, string-string and self-intersections result in formation of loops which subsequently decay. 
The energy loss rate can be approximated by [43] 

P--2^p-f, (2) 
a L 

where ' = d/dt and a{t) is the scale factor. The first term accounts for the expansion of the universe, and the second 
for string interactions with associated loop-formation. The network evolves towards a scaling regime, in which L is 
constant relative to the horizon dn ^ t. Indeed, setting L{t) = £,{t)t, it follows from Eq. (2) that 

where a{t) ^ (/3 = 1/2 in a radiation era, j3 = 2/3 in a matter era). The attractor scaling solution of Eq. (3) is 

C=[2(l-/3)]-i. (4) 

The "velocity- dependent one-scale model" (VOS) is slightly more sophisticated in that it introduces a dynamical 
velocity component v to the equations, describing the root mean square (rms) velocity of string segments. This model 
agrees quantitatively with results obtained from NG numerical simulations [25, 26]. The relevant equations are 

p = -2-il + v^)p-—^, (5) 
a L 

v^{l-v^)(^^-2%>j, (6) 

where the constant c represents the efficiency of loop formation, and k is the curvature parameter which indirectly 
encodes information about the small-scale structure on strings. It can be expressed as a function of the velocity [26]; 
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which incorporates the Virial condition v'^ < ^, observed in simulations. Following the same procedure as before, one 
now finds a scaling solution with 

2 _ k{k + £) 

^ 4/3(1 - /3) ' ^ ' 



2 



I3ik + £) 



(9) 



Of particular note for the following is that the velocity also enters a scaling regime (9) in which it stays constant in 
time. 

The relation between intercommutation probability V and the loop chopping efficiency parameter c is not fully 
understood at present. Nambu-Goto simulations of strings interacting with a microphysical probability V < 1 suggest 
that c ~ in both the matter and radiation era [45]. A different dependence, c ~ 7^^^^, was reported in [46] 

based on a flat space simulation. In the subsequent sections we follow [45] and take c to scale as the cubic root of the 
corresponding intercommutation probability. Such a weak dependence of c on P can be attributed to the presence 
of small scale structure on long strings, allowing for multiple chances of intercommutation when two string segments 
cross [45]. 



B. Evolution of multi-tension networks with junctions 



To describe the evolution of a multi-tcntion string network (MTSN) with junctions we adopt the model developed 
in [27, 36] in which we solve for the energy densities and rms velocities of each string type using the following set of 
equations: 



Pi 



-2^(1 
a 



E 

6, a<b 



LlLl 



(10) 



a 



E' 

6, a<b 



Vab {P.a+ Pb - Pz) ilb{t)Lf 



Pi 



r2 r2 
^a^b 



(11) 



Here /i^ is the tension of the ith type of string, and, as in the case of single tension networks, we define a correlation 
length Li through 



Pi 



Pi_ 



(12) 



In analogy to the single string case, the coefficients Ci in (10) quantify the efficiency with which self- interactions 
of strings of type i produce closed loops, thus removing energy from the long string network. The last two terms in 
Eq. (10) model the effect of collisions between strings of different types, leading to the formation of new segments 
ending on 3-string junctions. In particular, the penultimate term describes the loss of energy, from network i, due to 
string segments of type i colliding with segments of type a and forming links of type k. Similarly, the last term models 
the energy gain in network i through collisions between different strings a and b, leading to the formation of a link of 
type i. The parameter d^j = dj^, which we will discuss in more detail below, is essentially the probability with which 
strings of types i and j interact and produce a type k segment. This parameter captures interactions at the quantum 
level and incorporates volume effects [11, 21], as well as the kinematic constraints discussed in [28, 30, 32, 36]. The 
average length of the links formed by this process at time t is denoted by t^j{t), whose explicit form will be given 
below (Eq. (22)). In Eq. (11), the coefficients ki are curvature parameters which indirectly encode information about 
the small-scale structure on the strings. We will follow [26] and take them to depend on the rms velocities as 



2^2 /I - 



(13) 



^ab)^ 



The parameters b]^^ have been introduced in order to interpolate between the model of [27] (where 6^;, 
which the energy liberated by the formation of junctions is redistributed in the network as kinetic energy, and a 
model analogous to that in [24] (corresponding to b\f^ = 0), in which all of this energy is assumed to be radiated 
away. It is more likely that realistic situations are somewhere inbetween, with < c?^;,, corresponding to a fraction 
of the liberated energy being redistributed, the rest being radiated away. Finally, Vab is the magnitude of the relative 
velocity between strings of type a and b averaged over all directions, that is Vab = •\/ -t- . 
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C. Parameters for cosmic superstring networks 

In cosmic superstring networks, the tension of each string type is determined by the corresponding charges, (p, g), 
and the string couphng gg. In flat spacctime [17, 41, 42], 

M» = = ^Jphs + if : (14) 

9s 

where fip the tension of the lightest fundamental string (F string) carrying charge (1,0). The D string has a charge 
(0, 1), while the strings carrying charges {p,q) with p,q > 1 can be thought of as bound states between p F strings 
and q D strings. There is an infinite hierarchy of such {p,q) bound states, but, as was found in [24, 27, 36], the 
cosmological evolution of interacting networks of this type lead to solutions in which only the first few lightest strings 
dominate, allowing us to truncate the system at finite N . For our simulations we take N = 7 implying that we have 
seven different types of string carrying charges {pi, qi) = {p, q)i with 

{[p.q)^} = {(1,0), (0, 1), (1, 1), (2, 1), (1, 2), (3, 1), (1, 3)} , = 1, 7) . (15) 

The key parameters in Eqs. (10-11) are the self- and cross- interaction coefficients, q and d^^ respectively. They 
are averaged network quantities that depend on the microphysical intercommuting probabilities of the corresponding 
interaction processes, which in turn can be modelled using string theory techniques [21. 23]. The cross-interaction 
coefficients d^^ also depend on kinematic constraints on 3-string junctions [28, 30, 32]. 

In the remainder of this subsection we discuss in detail the dependence of the Ci and d'^^ , as well as other parameters 
entering Eqs. (10) and (11), on the fundamental parameters of string theory such as g^. Solutions of Eqs. (10) and 
(11) together with a discussion of the resulting network evolution is given in section II D: a reader less interested in 
the following more technical discussion may turn directly to section II D. 

1. Intercommuting probability Vij 

We begin by discussing the microphysical intercommuting probability Vij arising from interactions between strings 
of type i and j. From Vij one can determine both the self- and cross- interaction coefficients, c; and as discussed 
below. 

For processes involving at least one F string, that is F-{p, q) interactions, this probability Vij can be calculated 
perturbativcly in string theory [21, 23]. The result is a function of the string coupling g^, the relative velocity v and 
angle 6 of the incoming strings. However, there is also a model-dependent volume factor which depends on the size 
of the compact extra dimensions, parameterised by w (which will be defined precisely later), and by the amplitude 
of fluctuations of the string position fields. The latter depends on the string tension and, therefore, on the string 
coupling gg. A convenient way to write the answer is 

V^J{v,e,w,gg) = T.ij{v,e,gs)Vtj{w,gs) , (16) 

where J-ij{v, 9, g^) is the volume- independent quantum interaction piece and Vy(io, gs) encodes the volume dependence 
[21, 23]. 

la). The volume-independent piece, Tij{v,6,gs) 

As mentioned above, for interactions involving at least one F string, Tij{v,9, gs) corresponds to a perturbative 
amplitude that can be readily computed in string theory [21]. On the other hand, for interactions involving only D 
strings or heavier composites, the process is non-perturbative and less well understood: at present there are at least 
two approximate results, by Jackson et al. [21] (hereafter JJP) based on a worldsheet calculation, and by Hanany & 
Hashimoto [22] (hereafter HH) using a field theory approach. The two calculations are in good qualitative agreement, 
but there are quantitative differences reflecting the uncertainties arising from the currently incomplete understanding 
of such non-perturbative processes. Nevertheless, these calculations provide a basis for systematically computing the 
network coefficients in Eq. (10-11), and allow us to study the effect of these uncertainties on the macroscopic evolution 
of the networks (section II D). As the available methods for calculating these processes improve (for recent progress see 
[23]) we anticipate that such uncertainties will be eventually controlled. Here, we shall use the results of JJP for the 
perturbative processes involving at least one F string, and the result of HH for the non-perturbative D-D interactions. 
These have non-trivial dependencies on v and 6 and the relevant string charges p and q (for details see [21, 22]), but 
there are some key features with respect to their dependence on the string coupling gs which can be summarised: 
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• F-F string interactions scale with g'^, 

• F-(p, q) interactions with q > 1 scale with g^, 

• the non-perturbative D-D interactions scale with U^/ss ^ YfYieie is a number of order unity. 

For interactions between heavier composites, i.e. {p, q)-{p' , q') with q,q' > 1 and p,p' > 1, the amplitude is not known 
but it is understood that it is enhanced with respect to the D-D amplitude by the multiplicity qq' of the relevant Chan- 
Paton states [21]. In addition, for small values of the coupling gs we can neglect the effect of the light perturbative F 
strings so we will approximate the enhanced amplitude as J^{p.q){p' .q') = 1 ^ (1 ^ J^dd)'''' ■ The detailed form of the 
factors J-ij{v,9, gs) assumed in this paper (including their full v, 9 dependence) is shown in table I (left column). 



Interaction {ij) 




Vij 


F-F 


2 (l-cos By/ 1-1)2)' 


w 


F-D 


li2+(cos Oy/l-v'^) 


rmn{wgs '^^'^ , Sw, 1) 


F-(p,g), q>l 


9^i'^+(9sP-cose^(l-i,2){9fp2+,2))- 
^" 8sinSu^(l-ti2)(g2p2_,_g2) 


min (p^ + g^.gj^)^''''™ , Sw, 1 


D-D 






mm{wg^'\ 1) 


(p,g)-(p',g') , q,q'>l 


1- (I-.Fdd)"'' 


min { [(/ + q'gr'Kp" + q"g7')f' 



TABLE I: The coefficients J-^j and Vij for different string interactions. 



lb). The volume factors Vij{w,gs) 

These volume factors (see Eq. (16)) arise because the strings are moving in a higher-dimensional space. Thus they 
can miss each other as they cross, leading to an overall suppression on the interaction amplitude that scales with the 
inverse of the volume of the extra dimensions. However, it has been argued [21] that this suppression effect may not 
be as important as originally anticipated [11], because the string position fields are worldsheet scalars - not protected 
by any symmetry - and should therefore be stabilised at a minimum of a potential well (see, however, [47, 48]) rather 
than exploring the compact orthogonal dimensions. In this case, there is still a volume effect arising from the fact that 
strings are quantum objects whose positions fluctuate around the classical minimum, thus giving rise to an effective 
volume that each string explores, which is, however, a small fraction of the total. The size of the fluctuations is 
determined by the mass of the string, so this effective volume depends on the type of string and on the coupling gs 
(for details see [21]). For F-F interactions, both strings have the same tension and fluctuate by the same amount, but 
for F-D interactions (and for small gs) the fluctuation of the heavier D string can be neglected, leading to a volume 
which is a factor of (\/2)^ = 8 times smaller^ than the corresponding F-F volume. Finally, for D-D interactions the 
volume is approximately a factor gs = {^/g7f smaller than the F-F volume [21]. The overall volume suppression factor 
becomes unity when the effective string volume becomes equal to the minimum "string-scale volume" Kiiin = (27r^Q;')'^. 
Thus, deflning our model-dependent parameter w as 

W = Knin/l^FF G (0, 1] , (17) 

we have for the relevant suppression factors Vff = w, Vfd = min(8iy, 1) and Vdd = min(w(7j"^, 1) as long as 
gs <C 1. Here, we generalise these volume factors for (p, g)-strings and for ffs ^ 1 hi a phenomenological way so as 
to reproduce the above limits. In particular, we assign a factor of (p^ -I- q'^gj^)^^'^ to each (p, (7)-string. We show the 
resulting Vij{w,gs) in table I (right column). This is simply a choice we make in order to be able to systematically 
calculate these suppressions for different string types, but is one which successfully reproduces the results of [21] in 
the appropriate limits. The dependence of our network scaling results on this choice will be discussed in section II D 
below. We are thus left with a residual mo del- dependent variable w that we treat as an external tunable parameter 
like gs- Note that the choice w ~ 1 for this parameter (corresponding geometrically to a compactification very close to 



^ Note that 6 is the number of extra dimensions. 
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the string scale) makes Vy ~ 1 for all strings, so the dependence of the intercommuting probability Pij on the string 
type and gs is determined only by the quantum interaction Tij{v,9,gs) in this case. 

For the macroscopic string networks we are interested in, string collisions happen perpetually as the network evolves, 
with a range of relative velocities and angles. We can therefore average out the velocity and angle dependence of the 
microphysical probabilities Vij{v,9,w,gs) = Tij{v,9, gs)Vij{w, gs), simply by integrating over a Gaussian velocity 
distribution - peaked at the relative scaling velocities^ - and a flat distribution on the collision angles. This yields 
the probabilities Vij{w,gs) which now depend only on our two free parameters w and gs- 

The next step is to turn these probabilities into the network coefficients Ci and d'^j appearing in the macroscopic 
evolution Eq. (10-11). 

2. Self -interaction coefficients Ci 

For self-interactions, numerical simulations of Nambu-Goto string networks with reduced microphysical probabilities 
[45], suggest that the effective (loop-chopping efficiency) coefficient c in Eq. (5) scales with the third root of the 
microphysical probability V (refer to the discussion in section II A). Thus we will take here 

c. = S X Vl'^ , (18) 

where we have denoted Vu = Vi and c is chosen so as to reproduce single network results for Vi = \, d^^ = 0. Following 
[25, 26], we use c = Cr = 0.23 in the radiation era, and c = Cm = 0.18 in the matter era. 

3. Cross-mteraction coefficients d^j 

For the cross-interactions (i.e. with i ^ j) producing zipped configurations, there are at present no network sim- 
ulations to compare to in order to determine the dependence of d^j on the microscopic probabilities Vij, although 
for recent progress towards this direction see [33]. However, we may expect this dependence to be similar to the 
self-interaction case if these cross-interactions are initiated at a point (at which the incoming strings first cross) and 
then proceed by the zipping of the colliding strings, rather than exchange of partners. In that case, we may expect 
the effect of small-scale-structure on the strings to enhance the interaction probability by increasing the number of 
encounters within one crossing time, much like in the case of self- interactions'^. We therefore define: 

4=d.,5^, (19) 

with 

d, j =KX Vl^^ , (20) 

where k is a constant of order unity (we will set it to 1). Note that the effect of a k < 1 will still be captured to 
some extent^ by our subsequent analysis, because this parameter is degenerate with w (refer to table I and Eq. (16)), 
which we will vary as an external parameter. The scaling (20) is another model-dependent choice. Its effect will be 
examined in section II D. 

The final ingredient that enters the systematic computation of the cross-interaction coefficients d^j is the factor 

S^j in Eq. (19). This describes the conditional probability that the crossing of strings i and j produces a zipper of 
type fc, given that strings i and j have interacted. The interaction is fully described by an additive and a subtractive 
channel, so, for each pair {i, j} = {(p^, g^), (pj, qj)} there are only two possibilities for fc; either k ~ (pi + pj, qi + qj) 
or fc = {pi — Pj , qi — qj). Which channel is followed is determined by energetic considerations based on the balance of 
string tensions at the 3-string junction. These are nothing other than the kinematic constraints studied in detail in 
references [28-30, 32], and which must be satisfied for the junction to form in the first place. Reference [36] showed 



^ Note that this introduces an implicit dependence of the coefficients of Eqs. (10-11) on the actual scaling solutions for the string velocities 

Vi. This can be treated iteratively, by choosing initially some guess values for the scaling velocities to estimate these coefficients, and 

then update them iteratively using the returned scaling velocities until convergence. 
^ In an alternative picture, in which zipping would require the alignment of the incoming strings along a significant part of their correlation 

length, we would likely have a different dependence on small-scale-structure than in these point-like interactions. 

Since w G (0, 1] it cannot account for k > 1. 
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how these microphysical constraints can be integrated over a distribution of velocities and angles in a string network 
to obtain the averaged network coefhcients^ : 

4 = 1/ 9)) exp[(z; - v^,f/a.^] sin(0) v^dOdv , (21) 

where Q{—fp{v,9)) is a step function imposing the kinematic constraints fp{v,9) < (given by Eq. (51) of [30]), cr,^ 
is the variance of the velocity distribution, which is assumed to be a Gaussian peaked on the relative scaling velocities 
Vij = (vf + v^Y^'^ ^ and 5 is a normalisation factor (for details see Eq. (16) of [36]). Note the direct dependence of 
these coefficients on the string tensions /i^, ^j, /ife, and thus on the string coupling gg. By comparison, the dependence 
on w is very weak and enters only indirectly through the scaling velocities in Vij . 

In table II we show these suppression coefficients for the three lightest string components (F=l, D=2, FD=3) and 
for different values of the string coupling gs in the radiation era^. Due to the weak dependence of these coefficients 
on w, there is no significant difference between their values at it; = 0.1 and w = 1. 



w 




c3 


c2 
■""IS 


S23 






= 0.04 


0.180 


0.293 


0.950 




gs 


= 0.1 


0.117 


0.302 


0.881 




9s 


= 0.2 


0.071 


0.312 


0.790 


w € [0.1, 1] 


9s 


= 0.3 


0.050 


0.325 


0.707 




9s 


= 0.5 


0.033 


0.354 


0.590 




9s 


= 0.7 


0.028 


0.388 


0.516 




9s 


= 0.9 


0.026 


0.424 


0.462 



TABLE 11: The coefficients S^j for the three lightest string components (F= 1, D= 2, FD= 3) for different values of the free 
parameters ga and 0.1 < w < 1, in the radiation era. 

Putting all these factors together, the resulting coefficients Ci and df^-, in the radiation era, for interactions between 
the lightest string components are shown in table III for the same values of gs and w as used in table II. In the matter 
era, the coefficients Ci are a factor 0.78 smaller, while the d^j coefficients change only through the implicit velocity 
dependence in Eq. (21). This change in d^^ does not affect significantly the scaling values of the string correlation 
lengths and velocities, and so can be neglected. The main difference with reference [36] is that in [36] the coefficients 
dij were chosen independently as free parameters, while now they are systematically computed as described above 
and they depend only on and w, which are our free parameters. 



w 


9s 


Cl 


C2 


C3 


d?2 








w 


gs 


Cl 


C2 


C3 


rf?2 


rf?3 


rf23 




9s 


= 0.04 


0.02 


0.13 


0.13 


0.05 


0.08 


0.55 






gs 


= 0.04 


0.01 


0.13 


0.13 


0.05 


0.07 


0.55 




gs 


= 0.1 


0.03 


0.16 


0.16 


0.04 


0.11 


0.62 






9s 


= 0.1 


0.02 


0.16 


0.16 


0.04 


0.10 


0.62 




gs 


= 0.2 


0.05 


0.19 


0.19 


0.03 


0.14 


0.63 






9s 


= 0.2 


0.02 


0.19 


0.19 


0.03 


0.13 


0.63 


w = 1 


gs 


= 0.3 


0.07 


0.20 


0.20 


0.03 


0.16 


0.61 




w = 0.1 


9s 


= 0.3 


0.03 


0.20 


0.20 


0.02 


0.14 


0.61 




gs 


= 0.5 


0.10 


0.21 


0.21 


0.02 


0.21 


0.54 






9s 


= 0.5 


0.05 


0.20 


0.21 


0.01 


0.15 


0.54 




gs 


= 0.7 


0.12 


0.22 


0.22 


0.02 


0.26 


0.49 






9s 


= 0.7 


0.06 


0.15 


0.22 


0.01 


0.17 


0.39 




gs 


= 0.9 


0.15 


0.22 


0.22 


0.02 


0.31 


0.45 






9s 


= 0.9 


0.07 


0.12 


0.21 


0.01 


0.20 


0.31 



TABLE III: The coefficients Ci (in radiation era) and djj of Eq. (10-11) for the three hghtest string components (F=l, D=2, 
FD=3) for different values of the free parameters Qs and w. The values of a in the matter era are obtained by multiplying the 
Ci in this Table by a factor Cm/cr = 0.78 (see Eq. (18)). Finally, the values of d*j in the matter era change only through the 
implicit velocity dependence in Eq. (21) and this change does not affect significantly the scaling values of the string correlation 
lengths and velocities. 



Notice that we have corrected an error in reference [36], where the kinematic constraints were effectively overcounted due to the 
appearance of an extra probability P^*, see Eqs. (4) and (18) in [36]. 

There is a wealc dependence of the coefficients S^- on the expansion law (matter vs radiation era) due to their implicit dependence on 
the scaling string velocities through the Gaussian distribution in Eq. (21). 
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4- Average length of zippers 



The average length £ij{t) of zippers produced at time t which appears in Eq. (10-11) is, in principle, model- 
dependent. However, for cosmic superstrings, which are not subject to topological conditions which exist for example 
in non-Abclian field theory strings, it can be taken to be [27, 36] 

With this choice we are assuming that the produced zipper has a length which is smaller than - but close to - the 
smallest of the two correlation lengths of the colliding strings. 



5. The bij coefficients in the velocity Eq. (11) 

For simplicity we set the coefficients 6jj to zero. In this way we concentrate on the effects coming from the junction 
terms in Eq. (10) and not from the extra terms in the velocity evolution Eq. (11). In section II D we will also discuss 
the case with = d^j (refer also to section II B). In agreement with [36] we find that our key results are insensitive 
to this choice. 



6. Radiation to matter era transition 



Finally, when solving Eqs. (10-11) we must of course run through the regimes of radiation domination to matter 
domination and finally to a A dominated era. Of particular importance is the transition from radiation to matter 
domination and we interpolate between these two eras by following the approach in [37, 38], namely by replacing, for 
each string type, 

= X rl'\ (23) 

Here Cr = 0.23 (radiation), c,„ ~ 0.18 (matter), g = 300, and the scale factor a{t), calculated numerically from the 
Friedmann equation, is normalised so that a = 1 today. 



D. Scaling of cosmic superstring networks at large and small string couplings 

We have solved Eq. (10-11) numerically using the parameters computed in section II C for different values of our 
external parameters gs and w. As in references [24, 27, 36], we find scaling solutions in which all network components 
reach approximately constant string number densities and rms velocities during the radiation and matter eras. 

The string number density per unit Hubble volume is defined by 

^ {t/L,)\ (24) 

We find that it is dominated by the lightest three network components, namely the F, D and FD strings, while all 
heavier components generally have negligible number densities. The same three strings also dominate the effect on 
the CMB obscrvables or, more generally, any observable linearly related to the two point function of the string energy 
momentum tensor. With that in mind, it is useful to introduce the power spectrum density Mi given by [4] 

M.^d)". (25) 

which determines the amplitude of string induced power spectra. 

Fig. 1 shows the evolution of the rms velocity, number density and the power spectrum density for the three lightest 
strings for gs = 0.04 and gs = 0.9, with w = 1. When gs is close to unity, the tensions of the F and D strings are 
comparable. As a result, their densities are comparable too, as shown in the panels on the right. When gs is decreased, 
the lighter F strings become more populous and thus dominate the number density of the network. However, the D 
strings become heavier at smaller gs and, while less numerous, can dominate the power spectrum density, as evident 
from the gs = 0.04 case shown in the left panels of Fig. 1. This transition from the power spectrum being dominated 
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FIG. 1: Evolution of the rms velocity Vi (top panels), number density = (t/Li)^ (middle panels), and the power spectrum 
density (fJ-i/^i)'^ (bottom panels) of the three lightest network components: F strings (solid black), D strings (blue dash), and 
FD strings (red dot), at two representative values of the string coupling gs. The panels in the left column are for Qs = 0.04, 
those on the right are for Qs = 0.9. All plots are for w — 1. For Qs — >■ 1, the tensions of the F and D strings are comparable, as 
well as their densities. At smaller Qs, the lighter F strings dominate the number density, while the heavier and less numerous D 
strings dominate the power spectrum. The epochs of radiation-matter equality and last scattering are indicated with vertical 
lines. Of particular relevance for CMB is the correlation length of the string type that dominates the power spectrum density 
at the time of LS (see also Figs. 2 and 3). The deviation from scaling at late times is simply due to the fact we are entering a 
A dominated era. 



by the light populous F strings at higher string couplings to heavy rare D strings at smaller couplings appears to be 
a generic property of cosmic superstring networks and is one of the key results of this paper. Reducing the volume 
parameter w to 0.1 results in an enhancement of the overall netvifork number density (since Vij oc Vij oc w), but the 
dependence of the network scaling patterns on remains the same. 

In Fig. 2, we plot the number and power spectrum densities of the three relevant strings at the time of last scattering 
vs gs for two values of w. Reducing the value of the string coupling from g^ 1, initially reduces the contribution of 
D strings to Mtotai as they become more rare, but this trend quickly changes as we continue reducing g^ and the D 
strings become heavy fijy ~ gj^, eventually dominating Mtotai in the case of w = 1. For w = 0.1, the power density 
Mjj increases with decreasing g^, but does not quite catch up with Mp over the range of string couplings we have 
studied numerically. Note, however, that the number density of F becomes a few orders of magnitude larger that that 
of D, meaning that the correlation length of F strings is very small. If their correlation length is smaller than the 
horizon size at last scattering, the F strings do not contribute significantly to the CMB polarisation and, despite their 
lower power density, the B mode spectrum will in fact be dominated by D strings. This will be discussed more in the 
next section. 

Having sketched this interesting trend in the dependence of network scaling patterns on g^, let us see the way 
it emerges from the functional forms for the MTSN parameters adopted in Section II C. At large string couplings 
{gs — > 1), the hierarchy in the string tensions is not pronounced {^b = SJ^/^f Mf), and the F and D strings 
can be comparable both in number and power spectrum densities. In particular, for w ^ 1 and gs ^ 1, we have 
Vf ''^ glw ^ Vjy ~ 1, so, in effect, one would observe a single effective network of tension ^ Md, with a correlation 
length similar to that of ordinary cosmic strings, but with the extra property of frequently forming 3-string junctions. 
For small w, the situation would be similar, but with a smaller correlation length, and the appearance of 3-junctions 
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FIG. 2: Dependence of the string number (top panels) and power spectrum (bottom panels) densities at the time of last 
scattering on the value of the string coupling gs for = 1 (left) and w = 0.1 (right), for the three lightest network components; 
F strings (solid black), D strings (blue dash) and FD strings (red dot). 



would be more rare since now Vf ~ ggW < w ^ 1 is somewhat smaller than ~ 'wg~^ > w, so D strings would be 
somewhat more rare. It should be noted, however, that at « 1 the perturbative methods used to calculate Vff are 
not expected to be accurate. 

For small gs, the difference between the tensions of F and D strings is large, and so is the difference in the 
corresponding coefficients Ci,d^^ (see table III). As a result, the much lighter F strings dominate the network number 
density. The power spectrum density, however, is dominated in this case by the less populous, but much heavier 
(/iD oc gj^/ip) D strings which, being very massive, evolve practically independently of the light F-string network. 
Indeed, even the zipping between F and D strings gives rise to FD composites with a tension practically equal to that 
of the heavy D strings. So, as far as the contribution to the power spectrum is concerned, one has effectively a single 
string network of D strings ~ the F string component remains unobservable due to its low tension, despite the fact 
it has a high number density. Moreover, unless w <^ g^, the volume factors Vdd wgj^ in Eq. (16) approach unity 
while the non-perturbative factor J^dd ^ 1- As a result, the network properties are similar to those of an ordinary'' 
field theory string network with c~V~ 1. For w <^ g'^ the situation is similar, but the number density is enhanced 
(smaller correlation length) compared to ordinary cosmic strings. 

There may also be an intermediate regime, in which the string coupling gg is small enough for the light F strings 
to dominate the number density, but large enough for the heavier D strings to contribute significantly to the power 
spectrum density without dominating it completely. In this case there is a notable difference in the coefficients cp and 
CD and there are two substantially different correlation lengths describing the network. Any possible observational 
effect from this hybrid network would depend on some combination of the two correlation lengths that could be 
significantly different than for ordinary cosmic strings. Such an intermediate regime, while potentially interesting, is 
not the main focus of this paper, and we will leave discussion of it for future work. 

In the next section we study possible observable signals induced by these networks on the CMBR temperature 



Note that for extremely small Qs it is possible that this behaviour changes again to one with cd ^ 1 as the non-perturbative amplitude 
is expected to fall to Tysd ^ 1 as — >■ (see [21]). This effect, however, is not understood quantitatively at present. 
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FIG. 3: The correlation length and the rms velocity at the time of last scattering (LS) as a function of the string coupling Qa 
for w = 1 (left) and w = 0.1 (right) for the three lightest network components: F strings (solid black), D strings (blue dash) 
and FD strings (red dot). We show the string type(s) that dominate(s) the power spectrum at LS with oversized points. The 
horizontal line at ^ = 0.05 indicates the thickness of the LS surface. 



and polarisation power spectra. Of particular relevance in this regard are the values of the correlation length (^ls = 
I/(tLs)/iLs) Eind the rms velocity (wls) of the string type that dominates the power spectrum density at the time 
of last scattering (LS). In Fig. 3 we show the dependence of ^ls and wls on gg- The fact that strings of different 
correlation lengths dominate the spectra at different values of gs means that there could well be distinctive signals in 
the CMB spectra as a function oi gg- This may allow us to use CMB spectra as a discriminant between these different 
types of limiting behaviour, potentially providing an invaluable tool for constraining the underlying string theory, in 
particular the string coupling g^ and, more indirectly, the scale of compactification. We also note that, as seen in 
Fig. 3, there is only a minor difference in the values of scaling velocities of the dominant string types at different 
couplings. 

Before ending this section let us briefly discuss the dependence of our results on the various model-dependent 
assumptions and approximations we have adopted. In agreement with [36] we find that our results are practically 
identical when 6fj = and = djj-. We have also considered deviations from the functional form of the volume 
factors Vij for (p, (7)-strings (third and fifth entries of Table I, right). Since the networks are dominated by the lightest 
strings, such deviations only affect the heavy, rare species, so our scaling results are robust at the percent level, except 
if (7s — 1- An interesting model-dependent extension of our string evolution model is to allow for spontaneous decay 
(unzipping) of heavy strings into lighter ones, by adding the appropriate new terms in Eq. (10). The effect of these 
terms is to transfer length from heavy to light strings, but, again, since our networks are dominated mostly by F and 
D strings, with heavier composites being suppressed, these terms only give rise to corrections at percentage level in 
most of the parameter space we consider^. The most important mo del- dependent assumptions are the sealing of dij 



When the number density of FD strings is comparable to that of D strings, then spontaneous decay leads to a reduction of the FD 
string density in favour of D and F strings, and this effect is important, in agreement with the findings of Ref. [33]. However, for cosmic 
superstrings, the combination of Vij and 5*^ (refer to Eqs. (16) and (19-21)) generally impUes that the FD string density is smaller than 
that of D strings (top panel of Fig. 2) and this effect is negligible. The effect is strongest for small Qs, when the densities of FD and D 
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with Vij, Eq. (20), and the adoption of a non-perturbativc amphtudc !Fij for interactions containing D strings in both 
of the coUiding segments (fourth and fifth entries in Table I, left). Even though the two limiting behaviours discussed 
above are qualitatively robust, the transition point between them, as well as the values of the correlation length of 
the dominant strings in each regime, depend quantitatively on both of these choices. This highlights the need for 
investigating numerical simulations of strings with junctions and for further developing non-perturbative techniques 
for D string interactions. 

Finally, we comment on the significance of using the one-scale assumption to model the evolution of the network. 
To test the role of this assumption in setting the hierarchy of densities of different string types we have incorporated a 
second scale in our model, following the approach of [27], i.e. adding an extra scale ^, common to all string species, and 
modifying appropriately the evolution equations for the correlation length and the velocity. We have concluded that 
the effect of the second scale is not large - the number densities of the strings do not significantly change and, most 
importantly, the hierarchy stays the same. Another potential concern, which is of relevance to the CMB predictions, 
is that the one-scale assumption equates the correlation length along the strings with the interstring distance. This 
is only a reasonable approximation for networks with sufficiently large intercommutation probabilities {Vi > 0.1). 
Hence, we do not consider values of w and gs at which intercommutation probabilities of the dominant strings are too 
small. 



III. CMB TEMPERATURE AND B-MODE SPECTRA FROM MULTI-TENSION STRINGS 

Having obtained the MTSN scaling solutions at different values of string couplings, we would like to examine their 
imprint on the CMB temperature and polarisation spectra. First, we describe how we calculate CMB spectra for 
MTSN, then turn our attention to the shapes of the spectra at large and small values of gs- 

A. Modeling CMB with CMBACT 

To evaluate the CMB temperature and polarisation spectra sourced by multi-tension string networks (MTSN) we 
modify the publicly available code CMBACT [38, 39] so as to allow for strings with multiple tensions whose scaling 
is modelled by Eqs. (10) and (11). 

In CMBACT, the string network is represented as a collection of uncorrelated string segments, an approximation 
proposed in [49] and adapted for calculation of CMB spectra in [1, 38, 40]. In the unconnected segment model 
(USM), straight segments of strings are produced at some early time and given random/ uncorrelated orientations and 
velocities. At later times, a certain fraction of the number of segments decays in such a way as to match the number 
density given by a scaling model. The initial positions and orientations of the segments are drawn from uniform 
distributions, and the direction of the velocity is taken to be uniformly distributed in the plane perpendicular to the 
string orientation (longitudinal velocities are neglected). 

In the default version of CMBACT, the key parameters of the segments — namely their length, rms velocity 
and number density — are modelled using the VOS Eq. (5) and (6). The USM does not explicitly follow the loop 
distribution, however the energy in the loops is effectively included as part of the covariant conservation of the 
energy momentum of the scaling network. On their own, the straight string segments with open ends violate energy 
conservation. To remedy this, CMBACT enforces energy conservation by calculating the components Tpo and T,;j (with 
i ^ j) of the energy momentum tensor and then using the covariant conservation equation V^T^^ = to calculate 
Toi and Tu. 

We emphasise that CMBACT is not a means for gaining new insight into the evolution of cosmic string networks. 
Instead, it is a tool for evaluating CMB spectra for given one-scale parameters, such as the correlation length and 
rms velocity. In [50] , it was shown that the CMB spectra obtained from field theoretical simulations of Abelian-Higgs 
(AH) strings [51] are reproduced by CMBACT when the one-scale parameters measured in the simulation are used as 
input. Also, the CMB spectra obtained from the NG simulations of [52, 53] were compared to those from CMBACT 
with the one-scale parameters measured in [52] , and good agreement was once again obtained. The default version of 
CMBACT uses the VOS model with parameters tuned to match NG simulations of [26, 54]. 



strings are closer, but in this case the (number) density of F strings is higher by at least two orders of magnitude. Thus, the effect of 
spontaneous unzipping in this case, is at most an order unity enhancement of the D string density and our results (which depend on 
the relative density and tension between F and D strings) remain robust. 
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The shapes of the string-induced CMB spectra are mainly determined by the large-scale properties of the string 
network, such as the correlation length and rms velocity. The overall normalisation of the spectrum has a simple 
dependence on the string tension fji and ^ = L/t [4, 7] given by Eq. (25). The CMB temperature spectrum (TT) 
receives a contribution from the LS surface, for which the relevant scale is ^ at the time of LS. In addition, TT 
receives roughly equal contributions at each subsequent epoch (which is the mechanism by which strings can produce 
a scale-invariant TT spectrum on large scales) and hence the value of ^ is approximately the value measured during 
matter domination. CMB polarisation, on the other hand, is sourced at the time of LS and thus the normalisation of 
the B-mode spectrum is given by ^ls • 

In this work we generalise CMBACT to include unconnected segments of N different types. The lengths and rms 
velocities of each type are determined from Eqs. (10) and (11). As in the single-tension case, the Tq^ and Ta components 
of the string stress-energy are determined from the covariant conservation equation, which now takes the form 

N 

v''5^r;, = o. (26) 

The overall amplitude of the CMB angular spectra is approximately determined by 

N W / \ 2 

Cf™«^(xMtotal=^il/,=^^) , (27) 
i=l ^=l ^ 

and the shapes of the spectra will be set by the correlation length and the rms velocity of the most dominant population 
of strings. 

We should note that CMBACT has only been tested against simulations of single tension strings with no junctions 
and with an intercommutation probability of unity. However, as discussed in Sec. II, the scaling solution of our FD 
networks tends to fall into two categories depending on whether the string coupling is large or small. Namely, the 
energy density of the network is dominated by light populous strings for large values of g^^ and rare heavy strings 
at much smaller gg. In either case, the junctions are rare, and the bulk of the anisotropy is seeded by a single type 
of strings. Moreover, the intercommutation probabilities of the dominant string species are 0.1 or larger implying at 
most a factor of 2 reduction in loop chopping efficiency This justifies the use of CMBACT for modelling the CMB 
spectra from FD networks at least in the two limiting cases of large and small string couplings. 



B. The CMB spectra for scaling FD strings 



Before proceeding to discuss the CMB spectra sourced by FD strings, let us first note that cosmic strings cannot 
contribute more than 10% of the total CMB temperature anisotropy [50, 55-61]. To comply with this bound, we will 
adjust the fundamental (F) string tension to be such that 

fs — ^Itrings / ^total = ^'l ' (28) 

where we follow conventions of [7] to define 

2000 

C^^ = ^(2^+l)Cr. (29) 

Fortunately though, even with a marginal contribution to the TT spectrum, strings can be a prominent source of 
B mode polarisation. This is because strings, unlike inflation, are actively sourcing vector mode perturbations of 
magnitude comparable to the scalar perturbations (see [2-4, 6-8] for work on the subject). 

The dependence of the TT and BB power spectra on the correlation lengths and the rms velocities was extensively 
studied in [7] and more recently in [50]. The overall amplitude of the spectrum is approximately given by Eq. (27). 
The correlation length and the rms velocity of the main string type set the dominant momentum modes in the strings 
stress-energy, which in turn determine the position of the main peak. Larger string correlation lengths will move the 
peaks in the TT and BB spectra to lower i. The dependence of the peak position on the rms velocity is non-monotonic. 
The positions of the TT and BB peaks move to higher multipoles (smaller scales) for low and moderate velocities, 
but move to larger scales (lower £) for higher velocities. Also, larger values of v somewhat decrease the amount of 
BB power relative to TT power. The current bound on the fraction of string sourced CMB temperature anisotropy is 
only weakly dependent on the detailed shape of the spectrum, which in turn depends on the different types of strings 
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FIG. 4: The normalised TT (left) and BB (right) power spectra for gs — 0.04 (solid) and gs = 0.9 (dash) for w = 1, normalised 
to give fs = 0.1. Note that the smaller string coupling leads to a discernible move in the peak of the BB spectra to smaller £. 




10000 



FIG. 5: The normalised TT power spectra for gs — 0.04 (left) and Qs = 0.9 (right), including the individual scalar (S), vector 
(V) and tensor (T) contributions. 



involved. For instance, the bound on global strings, which do not lead to a pronounced peak in the CMB spectrum, 
is also approximately 10% [58]. 

In Fig. 4 we show the TT and BB power spectra for two values of the string coupling gs (solid black line for 
gs = 0.04, and dashed line for gs = 0.9). In Fig. 5 we show the normalised total TT power spectra for gs = 0.04 
(left) and gs = 0.9 (right), including the individual scalar (S), vector (V) and tensor (T) contributions for each case. 
The scalar and vector contributions are of similar magnitude, as expected in the case of cosmic strings. The spectra 
are normalised to give fs = 0.1 as described earlier, which translates into G^p = 1.8 • 10~^ for gs ~ 0.04, and 
Guf = 2.1 • lO-'^ for gs = 0.9. 

Wc can interpret the CMB spectra at different values of gs in the context of the scaling behaviour of the different 
string types discussed in detail in Sec. II. In the case gs = 0.9, the F and D strings are dominant and most populous, 
having almost the same tension. Their correlation lengths and velocities are similar as well, with values close to those 
of ordinary strings. Hence, their contributions to the CMB spectra are comparable, whereas the contribution from 
the heavier FD string is not as important. On the other hand, in the gs — 0.04 case, the F string is again light and 
populous, but the D string is 25 times heavier than the F string, and so is the FD string. Thus, despite being very 
rare, the D and FD strings dominate the CMB spectra. In addition, there is a small but non-negligible contribution 
from the heavier (2, 1) string. 

The right panel in Fig. 4 clearly shows the impact of the changed hierarchy between the three kinds of strings. The 
heavy D strings which dominate Mi for gs = 0.04 have a larger correlation length, which translates into a BB peak 
at smaller £. This offers a tantalising possibility for probing for small string couplings in the CMB. 

In Fig. 6 we again show the B-type polarisation spectrum predicted by our string models for a 10% and a 1% 
contribution to the total TT. We compare this to the contributions from gravitational lensing of the adiabatic E- 
mode polarisation into B-modes. Of note, especially in the region of high i, is the existence of a possible detection 
window, where the cosmic string signal would manifest itself as an excess over the expected lensing contribution. 
The Planck satellite may be able to see the excess if strings contribute at a level currently tolerated by data [62], 
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FIG. 6: Left: The B-type polarisation spectra due to cosmic superstrings assuming a 10% contribution {fa = 0.1) are plotted 
with solid {qs — 0.04) and dashed {ge = 0.9) black lines. The expected spectra for E to B lensing (blue dot line) and from 
primordial gravitational waves assuming a tensor-to-scalar ratio of r = 0.1 (magenta-dot-dash line) are shown for comparison. 
Right: The magenta dot line is the lensing prediction, the black solid line is the sum of the string and lens-sourced B-mode power 
for ga — 0.9 for fa — 0.01. Strings manifest themselves via the systematic excess power at high-^ over the lensing prediction. 
The sum of strings and lensing contributions is also plotted for fa = 0.1 for ga = 0.9 (red dash) and ga — 0.04 (green dot-dash). 
By measuring the location of the main peak we can rule out either the small or the large values of ga . 




FIG. 7: The position of the peak of the BB spectrum as a function of the string coupling ga for w — 1 (left) and w — 0.1 (right). 

however, it is likely that Planck's TT, TE, and EE spectra wiW place tighter bounds on strings even without B-mode 
measurements. The science goals of the ground based experiments, such as the advanced stages of QUIET [63] and 
POLARBEAR [64], include accurate measurements of the BB spectrum from lensing. If these science goals are met, 
they should have the sensitivity to detect the excess due to strings at a level of /s ~ 10"'^. If the string contribution is 
sufficiently large, it may be possible to detect the main peak and thus rule out large or small values of gs based on the 
position of the peak. More studies are needed to determine the minimum value of fs for which a particular B-mode 
experiment will be able to detect the position of peak [9, 65]. In principle, it may be possible to de-lense the B-mode 
polarisation map, taking advantage of the fact that the B-mode due to lensing is a rotation of the E-mode, and hence 
E and B modes are highly correlated [66, 67]. Some preliminary forecasts of expected constraints on cosmic strings 
from de-lensed B-modes were reported in [5, 9] but more work is needed to understand the feasibility of measuring 
the position of the string induced peak from de-lensed polarisation maps [65] . 

In Fig. 7 we show the BB peak location as a function of gs, ranging from 0.04 to 0.9. For w = 1, shown on the left 
plot, we see that the position of the BB spectrum as a function of the string coupling gs is decreasing with decreasing 
gs- This is consistent with our earlier discussion for the behaviour of the power spectrum density, which starts being 
dominated by the F strings and, after a transition, ends up being dominated by the heavy rare D strings. Their 
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correlation length increases with decreasing g^, and the BB peak moves to smaller (.. 

In the right plot of Fig. 7, which considers the w = 0.1 case, we also see that the BB peak position decreases 
with a decreasing g^. However, here it happens for somewhat different physical reasons. Namely, as we see in Fig. 2, 
the power spectrum density Mi in the more populous F strings dominates through the whole range of gs we have 
considered. Yet it is again the D string contribution that dominates the B-mode spectra. The explanation for this lays 
on the smallness of the for the F strings. For ^ < 0.05, the string correlation length is smaller than the thickness 
of the LS surface, and thus most of the power in that type of string does not contribute to the B-mode. Only a small 
fraction of the total power in F strings contributes, that generated on larger scales, but it is much smaller than the 
contribution of the D strings. The important issue arising here is that when the ^ of the dominant string becomes so 
small, the one-scale approximation is no longer reliable. This is a well-known problem, and further improvements for 
accurately calculating the CMB contributions from such strings must be performed in the future. 



IV. COMBINED CONSTRAINTS ON ^jlf AND FROM CMB AND PULSAR TIMING 

As mentioned in the previous section, the amplitudes of CMB two-point correlations do not separately constrain the 
string tensions and their densities. Instead, they constrain a combination of and given by Eq. (27). In particular, 
they do not differentiate between dense networks of light strings and rare heavy strings. In the case of FD networks, 
the relative abundances of different types of strings are controlled by the string coupling gg, with the bound on the 
CMB normalisation leading to different values of hf for different values of gs- Thus, the requirement that strings 
contribute no more than 10% of the total CMB TT power can be translated into a joint constraint on and gs 
shown with a solid line in Fig. 8. As mentioned in the previous section, the dependence of this bound on the detailed 
position of the peak in the spectrum is relatively weak. The derived bounds on and g,, depend more directly on 
the normalisation of the spectra given by Eq. (27). 

The degeneracy between fii and ^i, or in the case of FD networks, between fip and gs, can be partially broken 
if other types of observations become available. One example would be a measurement of the position of the string 
induced bump in the BB spectrum. Based on the results in Sec. Ill, if a peak is found at ^pcak ~ 610 ± 50, that 
would strongly disfavour gs > 0.1, while £poak ~ 750 ± 50 would disfavour < 0.3. In the future, more sophisticated 
simulations of CMB from MTSN based on specific compactifications can, in principle, make more accurate predictions 
for the dependence of the BB peak position on gs- 

A second way to break the degeneracy between fip and gs is to combine bounds from the CMB with bounds on 
gravity waves (GW) emitted by strings, such as those coming from pulsar timing experiments [68] (for other recent 
constraints coming from pulsar timing and lensing sec [69]) and direct GW searches by LIGO [70]. The GW bounds 
constrain the energy density in strings approximately given by the combination /i/^^ for each type of string. The fact 
that the functional dependences of the GW and CMB bounds on /x^ and are different implies that by combining 
the two probes one can, in principle, reduce the degeneracy between fip and gs- 

To illustrate this point, we follow the procedure presented in [50], where the authors calculated the bounds on the 
cosmic string tension from pulsar timing [68] (note that they are stronger than the ones coming from LIGO). For a 
network of a single type of string, the formula is [71] 

^_,Gt, (l-{vl^)/^\ (l + 1.4.T)3/2_i 



where x = ac^/{TGfJL), a is the loop production size, and fim is the total matter density relative to the critical 
density. They use parameters measured from the Nambu simulations to give ^rad, ('f^rad) ^"^^ = 0.3, F = 60. 

We generalise this formula to include the three types of string that dominate the energy density of the FD network, 
namely the F, D and FD strings. Setting c = 1 we write 

. 1.17 X 10- > ' J '-^J^\ il + lAx.f'^-l ^ (3^) 




?rad,i" , 

where i = 1, 2, 3 correspond to the F, D and FD string respectively. Also, Xi = a/{T Gfii), and for simplicity we take 
the a and F parameters to be the same for all types of string. For a given value of gs, we use the extended VOS 
model of Sec. II to determine the values of With those in hand, given the bound ftgh^ < 2 x 10^^, which is the 

most reliable published limit [68], we can use the relations /i^i = fip/gs and upd ~ iJ'Fygs+^ to find the bound 
on Gfip. As in [50] we consider two limiting cases of <C 1 and 2> 1 
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FIG. 8: Left: Current bounds on hf and Qs from CMB (solid-black line) and pulsars (short-dash-green and dot-red line). Right: 
The forecasted bound on fip and gs from CMB based on future BB corresponding to 0.1% in strings (solid-black line) together 
with the bounds from pulsars for the case Xi ^ 1 with a = 0.001 (short-dash-green line). 



Ccise a: <C 1: In this case, eq. (31) becomes 

n,h' = I • 1.4 . 1.17 X 10-4 ^ G^. i'-^T^r^] , (32) 

SO it is independent of a (as in [50]). The corresponding joint constraint on /if and gs is shown with red dotted line 
in Fig. 8 

Case b: Xi ^ 1: In this case, eq. (31) becomes 

i=l \ / 

As in [50], substituting Xi ~ a/{TGfii) results in a 1/a dependence of fip- The corresponding joint constraint on fip 
and gs for a = 0.001 is shown with a short-dash green line in Fig. 8. Note that the chosen value of a gives a bound 
of G/i < 5 • 10"^ for a 'usual' cosmic string with S^rad = 0.13 and v — 0.65 [50]. 

In the limit Xi ^ 1, the bounds on /i^ are weaker than the ones coming from the CMB. On the other hand, for 
Xi ^ 1 the bounds can be much stronger, because of the 1/a dependence. However, the uncertainties in the upper 
bound on flgh^ in this limit are much more severe [50]. Note that the shapes of the two pulsar bound curves in Fig. 8 
for the two limits (x,; <^ 1 and Xi 3> 1) are different because the dependence on /i^ is different. 

Future probes of the B-mode signal will be able to constrain cosmic strings down to a 0.1% contribution to TT, 
which translates into an order of magnitude tighter bound on fi p. In the right panel of Fig. 8 we show the corresponding 
expected bounds from the CMB. 

In the context of specific brane inflation scenarios, the cosmic string based constraints on gs and hf would need to 
be considered in jointly with other predictions of brane inflation, such as the spectral index Ug, the tensor-to scalar 
ratio r and the tensor index nx, as well as possible departures from Gaussianity in the distribution of the primordial 
fluctuations. We expect that such a comprehensive approach can lead to non-trivial constraints on details of brane 
inflation models, and more generally, the fundamental parameters of string theory. 
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V. CONCLUSIONS 

Using the MTSN scaling model developed in [36] , we have studied the evolution of a network of cosmic superstrings 
with the aim of identifying characteristic trends in their scaling properties at different values of the fundamental string 
coupling Qg. One particularly interesting trend that we have discovered is in the so-called power spectrum density, 
which controls the amplitude of the two-point function of the string stress-energy. As 5^ — s- 1 we sec it is dominated 
by populous light F and D strings where as as decreases to smaller values rare heavy D strings come to dominate. 

Using the MTSN scaling model in CMBACT, we evaluated the contribution of the FD networks to the CMB 
temperature and polarisation spectra. We found that the difference between the scaling patterns at high and low 
values of gs are manifest as a shift in position of the peak in the B-mode spectrum. In the one-scale approximation 
adopted in this work, the correlation length is equal to the average inter-string distance, thus string networks of 
higher (lower) number density have smaller (larger) correlation lengths. The correlation lengths, along with the string 
velocities, determine the position of the peak. This points to the possibility of constraining gs with future CMB 
experiments, possibly providing us with the first opportunity to date to constrain this fundamental parameter of 
string theory with observations. 

Most observables, at least those that rely on averaged properties of string networks, constrain a combination of 
the string tension fj, and their number density, determined by Thus, in principle, it is quite hard to distinguish 
between the effect of many light strings vs a few heavy ones. Measuring a particular peak position in the B-mode 
spectrum would be one way to partially break this degeneracy, as we have discussed in this work. In addition, one 
can explore the fact that different observables constrain different combinations of ji and ^. For instance, while the 
amplitude of the CMB spectra is determined by (/i/^)^, the energy density of strings is proportional to /i/^^. In this 
paper, we have shown how this difference can be explored in the case of FD networks to partially break the degeneracy 
between the fundamental string tension jjLp and the coupling gs by combining the CMB constraints with those from 
bounds on gravity waves (GW) . 

The main trends we have identified in this paper are largely independent of many of the details of the underlying 
string theory model, as well as the assumptions that went into the CMB calculation and the predictions for GW. 
However, follow up studies are needed in several directions in order to make firmer quantitative predictions. For 
instance, we need to improve our understanding of the string interaction rates for different choices of gs and w, 
especially in the non-perturbativc regime. At very small string couplings, the density of the dominant species becomes 
so high that the one-scale approximation is almost guaranteed to break down. In such cases, a more sophisticated 
model is needed to properly describe the scaling of the network and its prediction for the CMB spectra. Whether it 
will be possible to measure a peak at high £ in the B-mode spectrum is another interesting question which is addressed 
in an upcoming publication [65]. It will depend strongly on the resolution and sensitivity of the experiments, as well 
as our ability to clean the contribution from weak lensing. The GW bounds on FD strings depend on the presence of 
cusps [72], and on the loop size distribution, which is not fully understood at present, and should be revisited in the 
light of additional GW signatures coming from y-junctions [73]. Also, while in all string models considered so far the 
B-mode from the ordinary strings is sourced predominately by vector modes, the tensor modes (i.e. large scale GW) 
were never properly worked out (since it requires accounting for the backreaction) . It is our hope that the potentially 
very exciting opportunity for testing fundamental theory based on the general trends identified in this work will serve 
as an additional motivation for pursuing the remaining open questions. 
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Appendix: Approximate analytical solutions 



A surprisingly significant amount of analytic progress can be made with just a few assumptions about the general 
forms of the scaling solutions. Let us first revisit the usual one-scale case. In Sec. II A we have seen that in the VOS 
model the scaling constant values of the rms velocity v^^ns and ^ are given by Eqs. (8) and (9) when k is given by 
Eq. (7). To obtain an estimate of their values, let us make an assumption that ~ 1/2, which is based on evidence 
from the numerical solutions. It suggests that a pcrturbative solution to the equations could be obtained by writing 



2 1 

" =2-' 



where e ^ 1. It then follows from Eq. (7) that k ~ Q\/2e/TT and hence from Eq. (8) 

^ 7r/3 

indicating that ^ = 0{e). In fact, from Eq. (8) and Eq. (7) we have 



(34) 



(35) 



^(1 - 8v^) 



where 



l-8w6 + C(l + 8w6) 



(36) 



(37) 



2V2' " /? 

This can be solved for e using Eq. (34), leading to solutions in the radiation (/? = 1/2) and matter {(i = 2/3) eras 



^ A ^ 

Crad = and e^at = Y "6 ■ 



(38) 



Note this implies that C = 0{f) during the radiation era, where as in the matter era C ~ C'(e^). To give a flavour of 
the values these take, if Cr — 0.23 (i.e. Cr ~ 0.26) in the radiation era and Cm — 0.18 (i.e. Cm ~ 0.2) in the matter 
era, then 



erad = Y = ^-^^ ^ ^'■"'^ - ^-^^ "^^^ '''^d - 0-41 (Wrad 0.64) 



whilst 



(39) 



^mat 




^ 0.18 ^ 



0.5 and wf^^j 



0.32 



(Wn 



0.57) 



(40) 



Let us now consider the MTSN model of Sec. II B with parameters of Sec. II C. The tension of a (p, g) string is 
given by Eq. (14) and, for simplicity, we consider the first three tensions for the F, D and the (1,1) string: 



Ml 



A'(i,o) — A'F 

M(o,i) = ^^Fg's 



= M(i,i) = ^^i-pyi + gs 



(41) 
(42) 

(43) 



Of course many other types of string could be considered, however in practice only these three string types will be 
used as they will be the most populous. For = 1, /zi = /i2 = Ma/V^, where as if ^ 1 then /ii ^ /X2 ~ Hs- 

Next, we consider the evolution Eqs. (10-11) in Sec. II B with the coefficients 6^^ set to zero. Then, as with the case 
of a single type of string, the velocity Vi of the ith string satisfies 



21 



with ki given by Eq. (13). Generalising the definition of e in Eq. (34) we write for the velocity of the ith string 

^'^l- (45) 
where Ci ^ 1. Eq. (13) then leads to ki ~ 6-\/2ej/7r, hence from Eq. (44) it follows that at scaling 

Ci^^ = 0{e) (46) 

providing the analogous solutions to Eq. (35). Using given by Eq. (22), the equation of motion Eq. (10) for the 
string density reduces to 

Of course, if the interaction term = 0, then this corresponds to the case of the single type of string. Now, since 

Vab^^vl+vi^l + 0{e), (48) 
it follows that by considering the order of each term in Eq. (47), we require 

= 0(e) radiation era (49) 

and 

djjfc = O(e^) matter era. (50) 
For three string types, Eq. (47) becomes in terms of e^, 

/I 3\ , ~ , , Tre? / A C B \ 



A 




C 




£2(61 + 


^2) 


h 

£3(61 + 


£3) 


B 




A 

h — 




£3(62 + 


^3) 


ei(ei + 


£2) 


C 




B 

h — 




ei(ei + 


^3) 


£2(62 + 


£3) 



where 



6^2 ( ^ - + ^4 ~C2{l~e,)^^( -^-^ + _ ^^—^ ) (52) 

6.3 ( ^ - + e4 _ (73(1 - .3) = ^ f -T7^ + -T-^ - ^3 ,,.,^, ) (53) 



A = d?2«12, S = 43^^23, C = dj^Vis ■ (54) 



To solve the system of Eqs. (51 - 53) we start by considering the results from the numerical simulations which were 
in the radiation dominated era (/3 = 1/2). For gs ^ 1, and for almost all values of w, we find that 

ci < C2 ~ C3 and A ~ C < B , (55) 

as seen, for example, in Table III. Then, as a first approximation we set A = C = 0, ci = 0, Vab ~ 1 and /3 = i. 
In that case Eqs. (51)-(53) reduce to 

3ei ~C^^ -^4 L 7^—^) (56) 

2 V £263(62 + £3) 
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As we have set C2 — C3 then C2 = implying that 62 = £3 with Eq. (57) and Eq. (58) reducing to 

£2 = £3 = I ^ — — 1 > £2 (59) 

where 62 corresponds to the solution Eq. (38) in which the strings are not interacting. The effect of B (with A = C = Q) 
is therefore to reduce the velocity of the D and (1, 1) string, but to increase their correlation length and thus decrease 
their number density. From Eq. (56) 

3£i-Ci=-^ (60) 

where €2 is given in Eq. (59). Since ci ^ C2 we expect ei ^ £2 hence to a first approximation the solution for the F 
string is that corresponding to a network of only one type of string, namely 

. = f . m 

implying that the velocity of this string will be very close to l/\/2 and 6 will again be small. In terms of 6 and the 
original Ci we have 

6 = V2ci , 6 = 6- a/2c2 + ^3 (62) 

The comparison with the simulations is very encouraging. For example for the case w = 0.1, the prediction for 6 in 
Eq (62) matches the simulations exactly for 0.04 < gg < 0.5. The approximation we have adopted has forced 6 = 6 
and the predictions match very well the numerical solutions for 6 for values of < 0.3, differing by no more than 
7%. It fails though to allow for the observed deviation between 6 and 6 although this could well be accommodated 
for by an iterative approach to the solution. 
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